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We investigate cosmological scenarios in the theory of gravity with the scalar field possessing 
a non-minimal kinetic coupling to the curvature. It is shown that the kinetic coupling provides 
an essentially new inflationary mechanism. Namely, at early cosmological times the domination of 
coupling terms in the field equations guarantees the quasi-De Sitter behavior of the scale factor: 
a(t) oc 6 Hk * with H K = l/y/9K, where re ~ 10 -74 sec 2 is the coupling parameter. The primary 
inflationary epoch driven by non-minimal kinetic coupling comes to the end at tf ~ 10 -35 sec. 
Later on, the matter terms are dominating, and the universe enters into the matter-dominated 
epoch which lasts approximately 0.5HQ 1 ~ 0.5 x 10 18 sec. Then, the cosmological term comes 
into play, and the universe enters into the secondary inflationary epoch with a(t) oc e HAt , where 
Ha = \/A/3. Note that the present value of the acceleration parameter q — da/ a 2 is estimated 
as go — 0.25, that is the universe is at the beginning of the epoch of accelerated expansion. Thus, 
the cosmological model non-minimal kinetic coupling represents the realistic cosmological scenario 
which successfully describes basic cosmological epochs and provide the natural mechanism of epoch 
change without any fine-tuned potential. 

PACS numbers: 98.80.-k,95.36.+x,04.50.Kd 



I. INTRODUCTION 

In recent decades the observational cosmology has been going through the period of the rapid growth. Precise 
measurements of the Cosmic Microwave Background (CMB) radiation 1 , systematic observations of nearby and 
distant Type la supernovae (SNe la) [2] , study of baryon acoustic oscillations [5] , mapping the large-scale structure 
of the Universe, microlensing observations, and many other remarkable results (see, for example, the review [1]) have 
essentially expanded our knowledge about the Universe. Amazing discoveries, such as the accelerating expansion of 
the Universe and the dark matter evidence, have set new serious challenges before theoretical cosmology faced the 
necessity of radical modification of the standard model having successfully been exploited for a long time. Now, any 
viable cosmological model has to be able to describe several qualitatively different epoches of the Universe evolution, 
including the primary inflation, the matter-dominated stage, and the present acceleration (or secondary inflation). 
Moreover, it should also explain a mechanism providing an epoch change. These challenges have prompted many 
speculations mostly based on phenomenological ideas which involve new dynamical sources of gravity that act as 
dark energy, and/or various modifications to general relativity. The spectrum of models, having been postulated 
and explored in recent years, is extremely wide and includes, in particular, Quintessence [5], K-essence [6], Ghost 
Condensates [7], Dvali-Gabadadze-Porrati gravity 0, f(R) gravity [5], and others (see, e.g. Refs [TOUT?] for detailed 
reviews of these and other models). 1 

It is worth noticing that the most of phenomenological models represents various modifications of scalar-tensor 
theories. For example, the quintessence is the ordinary scalar field with a fine-tuned potential; the K-essence is a 
scalar-tensor theory with a generalized kinetic term; the f(R) gravity could be re- written in the Einstein frame as the 
ordinary general relativity with some effective scalar field, and so on. A broad class of cosmological models includes 
scalars non-minimally coupled to gravity p~8l420j . Additionally, one can further extend scalar-tensor theories, allowing 
for non- minimal couplings between the derivatives of the scalar fields and the curvature [21j . and these scenarios 
reveal interesting cosmological and astrophysical behaviors [2"2Tl26| . 

In our recent works [241 125) we investigated cosmological scenarios with the non-minimal kinetic coupling between 
the scalar field and the curvature, examining both the quintessence and the phantom cases in zero and constant 
potentials. According to the parameter choices and without the need for matter, we have obtained the variety of 
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1 This plethora of models reflects a deep crisis of phenomenological approach in the modern theoretical cosmology. To date, there are no 
unique criteria to prefer one or another phenomenological model. 
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behaviors including a Big Bang, an expanding universe with no beginning, a cosmological turnaround, an eternally 
contracting universe, a Big Crunch, a Big Rip avoidance and a cosmological bounce. However, the most interesting 
and important feature we have found is that the non-minimal kinetic coupling provides an essentially new inflationary 
mechanism and naturally describe transitions between various cosmological phases without any fine-tuned potential. 
The inflation is driving by terms in the field equations responsible for the non-minimal kinetic coupling. At early 
times these terms are dominating, and the cosmological evolution has the quasi-de Sitter character a(t) ~ e Hnt , where 
H K = 1/-\/9k and k is a coupling parameter. Later on, in the course of the cosmological evolution the domination 
of K-terms is canceled, and this leads to a change of cosmological epochs. In particular, in the cosmological scenario 
with zero potential the primary inflationary stage with a(t) ~ e H,it driving by non-minimal kinetic coupling changes 
into a power-law expansion with a(t) ~ i 1 / 3 [53], whereas in the case of constant potential V(<f>) = A/8tt the universe 
generally transits from one quasi-de Sitter phase to another with a(t) ~ e^ A *, where = \/A/3 25J. 

In this paper we consider more realistic cosmological models including an ordinary matter and the cosmological 
constant in addition to the scalar field with non-minimal kinetic coupling. 



II. ACTION AND FIELD EQUATIONS 



In this section we present the cosmological paradigm with non-minimal kinetic coupling between a scalar field and 
the curvature. In order to describe the quintessence and the phantom field in a unified way we adopt the e-notation, 
that is the parameter e takes the value +1 for the canonical field and —1 for the phantom one. 



A. General fields equations with non-minimal derivative coupling 



Let us construct a gravitational theory of a scalar field <f> with non-minimal derivative couplings to the curvature. 
In general one could have various forms of such couplings. For instance, in the case of four derivatives one could have 
the terms Ki-R</> jAt ^ ,M , K 2 R f _ tl y(j>' tJ '4>'' u , n 3 R(j)D(j) 7 K i R fJi ^(j)(f>' ,lu , n^R.^cfxp^ and K e DR(j) 2 , where the coefficients k 1; . . . , k 6 
are coupling parameters with dimensions of length-squared. However, as it was discussed in [2TJ |22J [23J , using total 
divergencies and without loss of generality, one can keep only the first two terms, KiRifi^c/)'^ and K^R^v^^ 11 . 
Generally, field equations with non-minimal derivative coupling are of third order. However, in the specific case 
— 2/ci = k 2 = k the coupling term G /i „</>' M </>' l/ gives field equations of second order [2~4"j . 

Further we will consider cosmological scenarios in the theory of gravity with the action 

S = J d 4 x^g~ | ^ - [e 9tw + kG^] W - 2V(<f>) J + S m , (1) 

where S m is the action for ordinary matter (not including the scalar field), V(<p) is a scalar field potential, g^ is 
a metric, R is the scalar curvature, G^ u is the Einstein tensor, and k is the coupling parameter with dimension of 
(length) 2 . Varying the action with respect to g^ v and 4> gives the field equations, respectively: 

= 8^ [T#0 + T$> + kQ^] , (2a) 
[eg"" + kGHV,^ = -V+ t (2b) 

where = dV(<f>)/dtfi, is a stress-energy tensor of ordinary matter, and 

T$> = s[V ^V v ^>- \ 9llv {Vct>f}- gill/ V{<t>), (3) 

-V M V„^- \{V4>) 2 G^ (4) 

+g„u[- iv Q v^v Q v^0+i(n^ 2 



leads to the differential 



Due to Bianchi identity V M G M1 , = and the conservation law V M T^™' ) = 0, Eq. ( |2a 
consequence 

V»[T$+ K B tu ,]=0. (5) 

Substituting Eqs. ^ and Q into |5]), one can check straightforwardly that the differential consequence ^ is 
equivalent to pbl. In other words, Eq. (|2b|) is a differential consequence of Eq. (2a 
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B. Cosmological equations with non-minimal kinetic coupling 



Let us consider a spatially-flat FRW cosmological model with the metric 

ds 2 = -dt 2 +a 2 (t)dx 2 1 



(6) 



where dx 2 is the Euclidian metric, a(t) is the scale factor, and H(t) — a(t)/a(t) is the Hubble parameter. Denoting 
the present moment of time as t , we have ao = a(£o) an( l Ho = H(t a ). Note that, without loss of generality, we 
can set ao = 1; in this case a(t) is dimensionless. Supposing homogeneity and isotropy, we also get 4> = <fr(t) and 



T, 



(m) 



diag(p,p,p,p), where p = p(t) is the energy density and p = p(t) is the pressure of matter. 



Assume that the universe is filled with a perfect fluid consisting of several non-interacting components with the 
energy density pi and the pressure Pi. Each component obeys the conservation law VT^H 1 = 0, which in the metric 
^ takes the form of the continuity equation 

Pi + 3H( Pi + Pi ) = 0. (7) 
For a given barotropic equation of state pi = WiPi the continuity equation is easily integrated: 

Pi=ma~ 3( - 1+a,i) , (8) 
where pto is the energy density of the z-th component at t = to- The total energy density and pressure are found as 



22 Pi = z2 Ptoa 



3(l+iui) 



Wip i0 a 



-3(1+1 



The general field equations ^ written for the metric ^ yield 

3H 2 = 4tt0 2 (e - 9kH 2 ) + 8ttV(» + 8np 
2H + 3H 2 = -4tt0 2 
eQ> + 3H<fi) - 3k(H 



e + n (2H + 3H 2 + AH^- 1 
I + 2HH4> + 3H 3 (f>) = 



^■nV{(j)) - 8ttp, 



(9) 
(10) 



(11a) 
(lib) 
(11c) 



where a dot denotes derivatives with respect to time. Here it is worth noticing that only two of the last three equations 
are independent, because Eq. (11c) is a differential consequence of ( 11a I and (lib). Note also that Eq. (11c) can be 
represented in the more convenient form: 

[(e - 3kH 2 )4>] ' + 3H(e - 3kH 2 )<P = -V4,. (12) 



C. The case V(4>) = const 

Hereafter we will discuss the particular case V(</>) = const. In fact, in this case the scalar potential plays the role 
of a cosmological constant, and so we will use the notation 

V{4>) = A (13) 

supposing that the cosmological constant A is non-negative, i.e. A > 0. In the case V(cf>) = const, when V$ = 0, Eq. 
(12) can be easily integrated: 

- E < 2A >" 2 (14, 



where A is a constant of integration. 

Let us define the following dimensionless density parameters: 



^miO = ^4>0 — > ^AO — „ , (15) 

Per Per Per 
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where p cr — 3Hq/8tt is the critical density. Using these notations and taking into account Eqs. ^ and (14), we can 
rewrite the first-order equation (11a) in terms of density parameters: 



H 2 = Hi 



A0 2^ „3(l+».) + „6( £ _ SkH 2 ) 2 



At t = to this equation reduces to 



o +Y^o »0 O (e-9 K ff o 2 ) _ 



(16) 



(17) 



The latter represents a constraint for parameters f2 mi 0j ^o> ^ao- 



III. COSMOLOGICAL SCENARIOS WITH NON-MINIMAL KINETIC COUPLING 



For given parameters, Eq. (16) completely determines the scale factor a(t) and hence the whole cosmological 



evolution of the Universe. Hereafter, we will assume that e = 1, supposing the normal scalar field with positive 
kinetic energy. Also we will assume that the perfect fluid filling the Universe has only one non-relativistic component 
with zero pressure p = and the energy density p = poa 3 - 



Model with k = and A = 0. 



First we consider the simple model with k = (no non-minimal kinetic coupling) and A = (no cosmological 
constant). In this case Eq. (16) reads 



h z = m 



with the constraint 



a" 



a" 



(18) 



^00 — 1- 



(19) 



It is easy to see that H oc a 3 — > oo if a — !► 0, and H oc a 3 / 2 — > if a — > oo. Integrating Eq. ( 18 ) together with the 



constraint (19) gives 



a{t) = [ln m0 Hl(t - t ) 2 + 3Ho(t - t ) + l] 



1/3 



(20) 



where f2 



vT 



E [0,1]. Note that the scale factor a(t) becomes zero at the cosmological singularity f* = t — 2(1 — 
^mo)/3^mo^o < to- Near the singularity a(t) oc (t — t*) 1 / 3 . In the particular case O m0 = 1 one gets the 
well-known solution a(t) = [^H (t — t ) + 1] 2 ^ 3 , describing the universe filled only with the dust. The opposite case 
O m o = corresponds to the universe filled only with the scalar field; in this case one finds a(t) — [3H Q (t — t a ) + l] 1 / 3 . 
Note that the scalar field behaves effectively like the stiff matter with p — p. It is worth also noticing that, in fact, the 
scalar field could play the role of dark matter. Actually, first, it represents an additional gravitational source being 
equivalent to the perfect fluid (stiff matter). In particular, on cosmic scales the scalar field is slowing down the rate of 
the cosmological expansion. Second, it interacts only gravitationally with another matter, and so one can speculate 
that dark matter halos of galaxies are formed from the gravitating scalar field. 

In Fig. [l]we give plots of a(t) for different values of fl m o- It is seen that the less is Q m o and, in turn, the greater is 
fi^o = 1 ~ ^mOj the slower is the expansion rate. This feature explicitly demonstrates the role of the scalar field as a 
dark matter component. 

Also, it is useful to represent graphically other cosmological characteristics: the Hubble parameter H = a/a and 
the acceleration parameter q = da/ a 2 , which is positive for accelerated and negative for decelerated expansion. The 
Hubble parameter H is determined by Eq. (18) as a function of a, i.e. H(a). It is easy to check that q(a) can be 
found as follows: 



(21) 
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where H' = dH/da. Plots of H(a) and q(a) are shown in Fig. [T] Note that q is negative for all values of a, i.e. the 
universe expands with deceleration in the course of the whole cosmological evolution. 

Summarizing, one can conclude that the simple model considered above describes satisfactorily the matter- 
dominated stage of the universe evolution. Moreover, the scalar field component could be regarded as a candidate for 
dark matter. However, this model is unable to describe the inflation and the late-time acceleration of the universe. 



B. Model with k = and A > 0. 



Now we take into account a positive cosmological constant A > in addition to the ordinary matter and the 
minimally coupled scalar field. In this case Eq. (161 yields 



H 2 = Hi 



^AO 



with the constraint 



AO 



raO 



a' 



(22) 



(23) 



reduces to (18) 

,) 1/3 



Therefore, the model with A > is also singular at t = t*, and near 
In the opposite limit a — } oo terms with il m o and fi^o are negligibly small 
^o^ao = A/3. The corresponding asymptotical solution has the 
quasi-De Sitter form a(t) oc e Hht , where H\ = i/A/3. 



In the limit a ^ Eq. |22| 
the singularity one has a(t) oc (t — 
in comparison with f^AOj an d so one gets H 2 



A general solution of Eq. (22 1 can be found in quadratures: 



flo(*-*o)= / 

J On = l 



a =i a [Q 



AO 



da 



^ a- 6 ] 1/2 ' 



(24) 



A graphical representation of this solution is shown in Fig. [2] for different sets of the density parameters f^Ao, Omo, 
and f^o- It is seen that at early stages the cosmological evolution has a matter-dominating character, whereas at late 
stages, when the cosmological constant is dominating, the universe expands with an acceleration. Fig. [2] represents 
also graphs for the Hubble parameter H — a/a and the acceleration parameter q — da/a 2 as functions of a. The 
behavior of q demonstrates that the universe expands with deceleration during the matter-dominated epoch, when q 
is negative; then q becomes positive and the deceleration is changed into acceleration signifying the beginning of new 
cosmological epoch. At the moment of epoch change the acceleration parameter is equal to zero, i.e. q = 0. Note 
that this moment depends on a ratio between the density parameters. In particular, the greater is J7ao the earlier is 
the beginning of the accelerating phase (see Fig. [2|. 

So, we can conclude that the model with a cosmological constant (or a constant scalar field potential) describes 
satisfactorily the matter-dominated stage and the late-time acceleration of the universe. However, the model is still 
unsuitable for describing an initial inflationary phase. 



(6) 



(c) 



FIG. 2: Plots (o), (6), and (c) represent graphs for a(t), H(a), and q(a), respectively, in the case n = and A > 0. Curves, 
from red to blue, are given for fi^o = 0.23, Qao — 0; 0.07; 0.27; 0.47; 0.73, and Q m o = 1 — ^<po — ^ao- 



Model with k > and A > 



Now, let us discuss the role of non-minimal kinetic coupling in the cosmological evolution. Assuming k > in Eq. 



( 16 1, we find 



h z = m 



^A0 



fyo(l - 9kH 2 ) 
a 6 (l-3 K H 2 ) 2 



with the constraint 



mO 



00 



(1 



3nH 2 ) 2 



(25) 



(26) 



In the limit a — > oo terms with f2 TO Q and f^o become negligibly small, so that Eq. (25) takes the asymptotical form 
H 2 pa H 2 fl\o = A/3, and the corresponding asymptotical solution has the quasi-De Sitter form a(t) oc e HAt . Note that 
it does not depend on the coupling parameter n at all. The essential role of non-minimal kinetic coupling manifests 
itself only at early stages. Actually, in the limit a — > the Hubble parameter has the following asymptotical behavior 
(see the appendix): 



H = a/1/9k + 0{a 3 



(27) 



Therefore, the early-time (t — > — oo) cosmological evolution has the quasi-Dc-Sittcr (inflationary) behavior with 
a(t) (x e ff »*, where H K = y/l/9/t. 

A general solution of Eq. ( [25] ) can be found numerically. In Figs.[3]a,& we represent plots of a(t) given for different 
values of n. In addition, we give graphs for the Hubble parameter H = a/a and the acceleration parameter q = da/ a 2 
as functions of a (see Fig. [3]c,rf). Graphical solutions explicitly demonstrate that a cosmological scenario in the model 
with the non-minimal kinetic coupling, positive cosmological constant and matter has three qualitatively different 
stages. At early times — oo < t < t\ the Universe is exponentially expanding, so that a(t) oc e fl " ! , H oc H K , and 
q oc 1. It is necessary to stress that the initial inflationary stage is completely provided with the non-minimal kinetic 
coupling of the scalar field to the curvature; the matter and cosmological constant do not play any role in this epoch. 
Then, at t\ < t < t 2 , the Universe expansion takes a power-law character. This is a matter-dominated epoch with 
the negative acceleration parameter q (see Fig. |3|. At late times t > t% the positive cosmological constant becomes 
dominating and the Universe expansion again takes the inflationary character such that a(t) oc e HAf , H oc Ha, and 
q oc 1. 



D. Estimations 



To estimate parameters of the phenomenological model constructed above, we will use well-established cosmological 



facts. As is known, the initial inflationary epoch should be over at tf 



10 



-35 



sec, and it should last at least 60 
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FIG. 3: Plots (a) and (6) represent graphs for a(t) in different scales; plots (c) and (d) represent graphs for H(a) and q(a), 
respectively, in the case k > and A > 0. Curves, from red to blue, are given for 7 = 3kHq = 0.001; 0.005; 0.01, Q m o = 0.04 
and Q$o — 0.23; the value of Qao obeys the constraint (26 1. The dashed curve in (a) corresponds to k = 0, Qao = 0, f2 m o = 0, 
£^<ao = 1; dotted curves correspond to n — 0, f^Ao = 0.73, f2 m o — 0.04, Q$ = 0.23. 



Hubble times (e-folds) (see, for example, Ref. [27]). Since a(t) oc e Hrit during this epoch, we have H K tj ~ 60, or 



6 x 10 d 



,1/2 



Then, taking into account that H K = W1/9k, we find k 



1CT 74 sec 2 , or 



10 



-27 



cm, where 



is the corresponding non-minimal coupling length. 



10 18 sec 1 , hence we obtain 7 



The present value of the Hubble parameter is Hq ~ 70 (km/sec) /Mpc ~ 10~ 18 sec -1 , 
10 -109 . Note that the value of 7 is extremely small in comparison with unity, and so neglecting terms with 7 in the 
constraint (26 1 yields Qao + ^mO + ^00 = 1- It is reasonable to assume Oao = ^dE) ^4,0 = ^dm, ^mo = ^fa> where 
Ode ~ 0.73 is the dark energy density, FIdm — 0.23 is the dark matter density, and fif, ~ 0.04 is the baryon density. 
With these values in hand we can calculate the acceleration parameter q(a) at the present time to with a(irj) = fflo = 1- 
Substituting f^Ao = 0.73, Q^o = 0.23, tt m o — 0.04 into Eq. (21), we find q = q(a ) = 0.25. It is worth noticing that 



<7o is positive, that is the model with given density parameters describes the present cosmic acceleration. Stress also 
that go < 1) what means the beginning of the secondary accelerated epoch. Later on, the accelerated expansion will 
enter into the quasi-De Sitter phase with q ~ 1. 

With given J7ao = 0.73, il^o — 0.23, and J7 m0 = 0.04 we can also estimate the duration of matter-dominated phase, 
i.e. tMD- Supposing that a = at the beginning and a = 1 at the end of the phase, we find 



HotMD 



da 



a [Oao + ^moa 3 + Q<t>oa 6 ] 



-61V2 



0.5. 



Hence, tMD — 0.5H Q 



0.5 x 10 18 sec. 



IV. SUMMARY AND CONCLUSION 



In this paper we have explored a cosmological model with the scalar field possessing non-minimal kinetic coupling 
to the curvature given as kG Ati ,</>'' J '4> ,p ', where k is a coupling parameter with dimension of (length) 2 . Additionally, the 
model includes an ordinary matter in the form of a pressureless perfect fluid, and the positive cosmological constant 
A, which could be also considered as the constant scalar potential V((f>) = A/8tt. It is worth noticing that two 
dimensional parameters k and A determine, eventually, two characteristic scales H K = l/\/9n and H^ = U A/3. 

The considered model represents the number of interesting and important features. 

• The most important feature is that the non-minimal kinetic coupling provides an essentially new inflationary 
mechanism which does not need any fine-tuned potential. The essence of the mechanism is that at early 
cosmological times the coupling K-terms in the field equations are dominating and provide the quasi-Dc Sitter 
behavior of the scale factor: a(t) oc e^"' with H K = 1/\/9k. Assuming that the initial inflationary epoch should 
be over at tj ~ 10~ 35 sec, and it should last 60 Hubble times (e- folds), one can estimate the coupling parameter 
k ~ 10~ 74 sec 2 and the corresponding length l K = k 1 ! 2 ~ 10~ 27 cm. 

• The cosmological scenario consists of three basic epochs: the primary inflation, the matter-dominated phase, 
and the secondary inflation. The primary inflationary epoch driven by non-minimal kinetic coupling comes to 
the end at tj ~ 10 -35 sec. Later on, the matter terms in the field equations begin to be dominating, and the 
universe enters into the matter-dominated epoch which is characterized by decelerated expansion. This epoch 
lasts approximately 0.5HQ 1 ~ 0.5 x 10 18 sec, then the cosmological term in the field equations comes into play, 
and the universe enters into the secondary inflationary epoch with a(t) oc e ffA *, where H\ = ^A/3. Note that 
the acceleration parameter q — da /a 2 during the inflation is q ~ 1. The present value of q is estimated as 
<7o ~ 0.25. This means that at present the universe is at the beginning of the epoch of accelerated expansion. 

• The model provides a natural mechanism of epoch change. Mathematically, every cosmological epoch corre- 
sponds with an appropriate term dominating in field equations in a definite time interval, and a change of epochs 
occurs when one dominating term is replaced by another. In this connection, it is necessary to stress that now 
one needs no reheating mechanism to create matter after the primary inflation. Actually, the ordinary matter 
exists during the whole period of the inflation driven by non-minimal kinetic coupling, however its influence on 
the cosmological evolution is suppressed. 

Summarizing, we can conclude that the theory of gravity with the scalar field possessing the non-minimal kinetic 
coupling to the curvature yields realistic cosmological scenarios which successfully describe basic cosmological epochs 
and provide the mechanism of epoch change. 
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Appendix 



Here we will discuss in details the equation ( 25 ) which reads 

n m0 , n 0o (i - 9k# 2 ) 



Using new notations 



we can rewrite it as follows: 



H 2 = Hi 



a 3 a e {l-3 K H 2 y 



y = (H/H ) 2 1 x = a 3 , ! = ZkH 2 , 



wi U2 1 - 3jy 
y = u)o-\ h 



x x 2 (1 — jy) 2 



(28) 



(29) 



(30) 
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The latter is a cubic equation for y(x). Let us discuss asymptotical properties of y(x). In the limit x — > oo we can 
neglect two last terms in Eq. ( 30 1 and find 



y = U) + O(x 1 ). 

To obtain asymptotics in the limit x — > 0, we represent y(x) as follows: 



(31) 



k=l 



Substituting this series straightforwardly into Eq. (30) yields three different asymptotical solutions at x — >• 0: 

4cJi 



(i) = 



37 277^2 



4 3w 2 + 47wf — 97Woo;2 2 
243 X 



0(x% 



y (2,3) = 7^1 ± V7 2 ^i - 127^2 | 
272 



(32a) 
(32b) 



Note that the first solution y^ has a regular behavior near x = 0; namely, j/W — > 1/37 if a; — >■ 0. It is worth noticing 
that the asymptotical value y^(0) = 1/37 does not depend on the parameters = {ujq, ui, uj 2 } at all, and so the 
asymptotic j/ 1 * 1 ~ 1/37 at x — )• has the universal character for Eq. (30). Note also that the solution is real for 
any real values of a;,. On the contrary, the solutions 1/ 2 ' 3 ) are singular at x = 0; their asymptotics are j/ 2,3 ) esc 1/x at 
x — > 0. Moreover, ?/ 2 > 3 ) could be real or complex conjugate depending on a sign of the operand of the square root in 
Eq. p2bl). 



To obtain a general solution, we rewrite Eq. ( 30 ) in an equivalent form 



y 3 + P2y 2 +piy + p o = 0, 



where 



P2 

P = 
Po = 



(2 + 7W )a; + 7W1 
72; 

(1 + 2 r yu!o)x 2 + 2ju>ix + 37^2 
7 2 a; 2 

LOqX 2 + LUlX + CJ2 

7 2 a; 2 



(33) 

(34a) 
(34b) 
(34c) 



are real. A cubic function y 3 + p 2 y 2 + Piy + Po with real coefficients Pi has generally three roots given by Cardano's 
formulae (see, for example, [28]): 



2/ 



2/ 

(2,3) 



J32 _ C _ p\ - 3pi 
3 3 



3C 



^ 3 )(x; 7 ,^)^-^ + 



P2 



C(l±n/3) , (lTiVS)(p 2 2 -3 Pl ) 



6 



6C 



where 



|(Q + 2pi-9p2Pi + 27po), 



9p 2 Pi + 27p ) 



4(P| 



3pi 



(35a) 
(35b) 

(36a) 
(36b) 



Note that the above expressions are not yet unambiguously defined, and so they need further explanation. To do this, 
we will assume that the square root refers to the principal (positive) square root when its operand is non-negative, and 
the cube root is to be interpreted as the real one, so that the solution is real, and ?/ 2 ' 3 ) are complex conjugate. 
Otherwise, there is no real square root and we can arbitrarily choose one of the imaginary square roots. Then, for 
extracting the complex cube roots of the resulting complex expression, we will choose among three cube roots that 
which provides j/W to be real. 



Since the root j/W is real, we will consider it as a physical solution of the cubic equation (33). It is necessary to 
emphasize that y^ 1 has the regular asymptotic (32) in the limit x — > (one can check this fact straightforwardly). 
Now, reverting to the old notations, we find 



H 



H 2 y (1) (a; 3kHq , fl AOl O m0 , tt^o) 



(37) 
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To obtain an explicit expression for H, one should use relations (34), (35), (36). However, a final expression is too 
cumbersome, and so we skip this step. Instead, we will discuss an asymptotical behavior of H. Taking into account 
Eqs. pll), (|32a|), we can easily find 



H = + 0{a- 3 ) if a^oo, 

H = y/l/9K + 0(a 3 ) if a^O. 



(38) 
(39) 



Finally, integrating Eq. ( 37 ) yields 



da/o 



>a =i vu- v-;3/«-ffoj fi AO,^mO,^</>o)] 1/2 ' 
This relation could be used for numerical calculations of a(t). 



(40) 
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